The magnetized discharge with dust: negative and positive space charge modes 
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The structure of a discharge across a magnetic field in a dusty plasma is analysed. The dust 
macroparticles are negatively charged, but are unmagnetized because of their high mass. The 
electrons are highly magnetized, and the ions have intermediate magnetization. This results in 
different transport rates of the different species across the magnetic field. Depending on the size 
of the magnetic field, and the relative charge on the different species, the dust grains can be the 
dominant current carrier. The space charge clouds near the electrodes will then be determined by 
the relative mobility of the different species. The discharge can operate in one of two modes, a 
positive space charge (PSC) mode, characterized by a strong cathode fall, and a negative space 
charge (NSC) mode, characterized by a broad anode fall. Features unique to the dust particles can 
also play a role in the structure of the discharge, such as the variable equilibrium charge on the 
grains, dependent on the local potential and species temperatures, the effect of gravity on the grain 
dynamics, and the rate of charging of the grains. The dust grains can also form an ordered structure, 
the dust-plasma crystal. A fluid model of the different species is used to calculate the structure of 
the resulting discharge, incorporating the above effects. The transition from the PSC mode to the 
NSC mode as the magnetic fleld, pressure and dust properties are varied is demonstrated. 

PACS numbers: 52.80.Sm, 85.10.Jz 



I. INTRODUCTION 

A Particle-in-Cell/Monte-Carlo (PIC/MC) simulation 
(van der Straaten et al. 1994, 1997) and semi-analytic 
treatment (Cramer, 1997) of the radial structm-e of a 
low pressure DC cylindrical magnetron discharge has re- 
vealed a potential and electric field structure highly de- 
pendent on the pressure and magnetic field. As either 
the pressure was reduced or the magnetic field was in- 
creased, the steady state discharge was found to exhibit 
a transition from a positive space charge mode (PSC), 
characterised by a strong cathode fall as occurs in an 
unmagnetised glow discharge, to a negative space charge 
mode (NSC) characterised by a broad anode fall. The 
reason for the transition to the NSC mode is the strongly 
reduced (according to classical theory) transport of elec- 
trons across the magnetic field in a low pressure, strongly 
magnetised plasma. These two modes of the magnetised 
discharge have been discussed by Thornton and Penfold 
(1978). 

There is little conclusive experimental evidence for 
the NSC mode. The cathode fall is always observed in 
planar magnetron experiments, even at very low pres- 
sures (Rossnagel and Kaufmann 1986, Gu and Lieber- 
man 1988). Experiments by Yeom et al. (1989) with a 
cylindrical magnetron with pressures and magnetic fields 
similar to those considered in this paper and in the PIC 
simulation study (van der Straaten et al. 1997) showed 
a distinct cathode fall and no anode fall over the entire 
range of discharge parameters. Langmuir probe measure- 
ments in a cylindrical magnetron were also reported by 
van der Straaten et al. (1997) with the same discharge 
parameters as used in the simulation, but the results did 
not agree with the simulation results in that no anode fall 
was observed at low pressures and high magnetic fields. 
However, Hayakawa and Wasa (1965) reported the exis- 



tence of a stable discharge operating in what appeared to 
be the NSC mode. The discharge featured a broad anode 
fall for a magnetic field strength greater than 4kG, which 
is considerably higher than the field strength predicted 
by the fluid model and the simulations (« lOOkG) for 
the onset of the NSC mode. In order to explain the per- 
sistence of the cathode fall in the experimental results it 
would be necessary for the electron transport across the 
magnetic field, at low pressures and high magnetic field 
strengths, to be considerably higher than is predicted by 
the classical transport coefficients. It has been postu- 
lated (eg Sheridan and Goree 1989) that turbulence or 
nonlinear coherent modes induced by instabilities in the 
partially ionized plasma in crossed electric and magnetic 
fields (Simon 1963) may increase the diffusion and drift 
of electrons, thus increasing their effective transport co- 
efficients. 

Dust macroparticles in a discharge are negatively 
charged, but are unmagnetized because of their high 
mass. The electrons are highly magnetized, and the ions 
have intermediate magnetization. This results in differ- 
ent transport rates of the different species across the 
magnetic field. Depending on the size of the magnetic 
field, and the relative charge on the different species, the 
dust grains can be the dominant current carrier. The 
space charge clouds near the electrodes will then be de- 
termined by the relative mobility of the different species. 
The two modes of the discharge will then be affected 
by the charge on, and the current carried by the dust 
grains. Features unique to the dust particles can also 
play a role in the structure of the discharge, such as the 
variable equilibrium charge on the grains, dependent on 
the local potential and species temperatures, the effect of 
gravity on the grain dynamics, and the rate of charging 
of the grains. The dust grains can also form an ordered 
structure, the dust-plasma crystal. A fluid model of the 
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Figure 1. The magnetron discharge geometry'. 



FIG. 1: The discharge geometry. 



different species is used to calculate the structure of the 
resulting discharge, incorporating the above effects. The 
transition from the PSC mode to the NSC mode as the 
magnetic field, pressure and dust properties are varied is 
demonstrated. 



dust. Thus we use Poisson's equation in one dimension, 
i.e. 



dE 
dx 



(3) 



where p is the net charge density and eg is the permittiv- 
ity of free space. This equation may be rewritten as 



dE ^ ]_ 

dx £0 



VTe 



(4) 



where Vi is the ion drift velocity in the x-direction (neg- 
ative in this case), Vd is the dust grain drift velocity, and 
Ji, Je and Jd are the ion, electron and dust current den- 
sities in the x-direction. 
The dust current density is 



Jd = ndVdQd = ndQVdoQd 



(5) 



where we assume a constant flux of dust particles, with 
initial density and velocity Udo and Vdo- The dust charge 
varies in the discharge due to the varying local potential. 

Writing the total current density a,s J = Ji + + Jd, 
and assuming \vi\ <^ \vd\, we have 
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II. THE BASIC EQUATIONS 


dx 













(6) 



We consider a one-dimensional model of a magnetron 
discharge between two parallel plate electrodes, as shown 
in Figure 1. 

We assume a one-dimensional steady state distribution 
so that all quantities are a function only of the distance 
X from the cathode surface. The magnetic field B is 
uniform and directed parallel to the cathode surface. The 
ions and dust grains are assumed to be unaffected by the 
magnetic field, i.e. the ion and dust Larmor radii are 
assumed to be much larger than the distance d between 
the anode and the cathode. 

The electrons have two components of drift velocity, 
VTe transverse to the magnetic field in the x-direction, 
and v±e perpendicular to both the electric field E and to 
B, where 



VTe = E 



and 



V±e 



(1) 



where E is the x-component of E (negative in this case), 
Wc is the electron-cyclotron frequency and v is the colli- 
sion frequency of electrons with background gas atoms. 
The resultant drift of the electrons is at an angle 9 to the 
X-axis given by 



(2) 



The basic equations used are those of Davies and Evans 
(1980), modified to include the magnetic field and the 



where je is the fraction of the total current density due 
to electrons. A boundary condition that can be applied 
is that the electron current at the cathode is due solely 
to secondary emission of electrons caused by ion impact 
on the cathode. The secondary emission coefficient 7 = 
del di at the cathode is assumed known. 

The second basic equation we use is the electron charge 
conservation equation, or ionization avalanche equation. 
The electrons drift through the background neutral gas 
at the angle B to the x-axis and ionize the neutral gas 
molecules, and electron avalanches are formed. These 
avalanches are therefore also inclined at the angle B to 
the X-axis. If the coordinate along this direction is the 
normalized electron current density in this direction is 
jc,e and the electric field in this direction is i?^, and the 
ionization equation may be written 



(7) 



where a is Townsend's 1st ionization coefficient 
(Llewellyn- Jones 1966), 



= APexp(-C(P/|£;cl)') 



(8) 



where A and C are constants depending on the gas, P is 
the gas pressure and s = 1/2 for a monatomic gas. Since 
C,= xj cosB, j(^e = je/ COS and = £' cos 6*, [T] becomes 



3 



dje , . 

= aje 

ax 



AP 

cos 6 



exp (^-C{P/\E\ 



\l/2 



Je. (9) 



The only difference in equation ^ to the unmagnetized 
case is therefore the replacement of the pressure P by the 
"effective pressure" P/ cos9. 

The ion mobility is assumed unaffected by the mag- 
netic field, so we assume, as do Davies and Evans (1980), 
that 
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(10) 



where fc is a constant. This gives a good representation of 
the experimental ion drift (Ward 1962). However we note 
that this means that the ion and electron drift velocities 
have different E dependences, so r is not strictly inde- 
pendent of E as we have assumed so far. A dependence 
of r on £^ would prevent the application of the analysis 
used here, so we neglect it, noting that it could cause an 
error in our results at high magnetic fields. 

The charge of a (negatively charged) dust particle is 
determined by the current balance equation 
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III. RESULTS 



The equations have been solved for a number of cases, 
using the above prescription, to illustrate the effect on 
the discharge of increasing the magnetic field and the 
density of dust particles. The parameters used in the 
numerical examples are those of van der Straaten et al 
(1994), viz. d = 2.2cm, P = 5mTorr and 50mTorr, 
and Argon gas, for which A = 29.22cm-iTorr-\ C = 
26.6V^/^cm^^/^Torr^^^^, electron mobility for zero mag- 
netic field fie ^ 3 X lO^cm^TorrVs"^ and k = 8.25 x 
lO^cm'^/^Torr^/^V^^/^s"^. The corresponding electron 
collision frequency is i' = 6 x lO^Ps^^ where P is in 
Torr. The ratio of ion and electron mobilities in the un- 
magnetized gas is ro — 3.3 x 10~^ (Ward 1962). 

Figure 2 shows the electric field, ion velocity and space 
charge profiles for a pressure of 5mTorr. The cathode is 
at a; = and the anode is at x = 2.2cm. No magnetic 
field and no dust is present, and the potential drop — 
209V. There is the usual positive space charge region 
near the cathode. 

Figure 3 shows the electric field, ion velocity and space 
charge profiles. No magnetic field is present, but dust is 
present, with n^o — lO'^i and the potential drop = 209V. 
A negative space charge region forms near the anode, due 
to the dust. 

Figure 4 shows the electric field, ion velocity and space 
charge profiles. No magnetic field is present, but dust is 
present with a higher density, with n^o = 3 x 10'^, and 
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FIG. 2: The electric field, ion velocity and space charge pro- 
files for a pressure of 5mTorr. There is no magnetic field or 
dust. 



the potential drop = 209V. The negative space charge 
region is more pronounced. 

In figure 5 a magnetic field is present, with cos 6 — 
0.2, but no dust is present, and the potential drop — 
392V. The negative space charge region is now due to 
the magnetic field. 

In figure 6 a magnetic field is present, with cos 9 = 0.2, 
and dust is present with n^o = 10^, and the potential 
drop — 392V. The negative space charge region is more 
prominent, due to the dust. Dust charges negative near 
the NSC. 
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FIG. 3: The electric field, ion velocity, space charge and dust 
particle charge profiles. There is no magnetic field, but nao = 
10^ 



In figure 7 a stronger magnetic field is present, with 
COS0 = 0.5, and dust is present with n^o = 10'^, and the 
potential drop = 524V. The negative space charge region 
is wide. Dust charges positive over most of the discharge. 
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FIG. 4: There is no magnetic field, but rido = 3 x 10^. 



IV. DISCUSSION AND CONCLUSIONS 

A numerical solution of the electron, ion and dust fluid 
transport equations for a magnetized discharge has been 
developed, building on previous work for an unmagne- 
tized steady-state glow discharge. Understanding the 
transport of charged particles across the magnetic field is 
important for modelling the operation of magnetron de- 
vices used in plasma processing for industry. The effects 
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FIG. 5: A magnetic field is present, with cos Q = 0.2, but no 
dust is present. 

of charge on the dust particles on the transition from a 
positive space charge mode to a negative space charge 
mode as the magnetic field is increased or the pressure is 
reduced has been demonstrated. 

The presence of dust can create a negative space charge 
region near the anode, which enhances or mimics the ef- 
fect of a magnetic field. If however the field is so strong 
that the ions are magnetized (future work), the dust 
grains may carry most of the current, which will enhance 
the positive space charge region. 
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FIG. 6: A magnetic field is present, with cos 6 — 0.2, and 
nao = 10^ 
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FIG. 7: A magnetic field is present, with cos 9 — 0.5, and dust 
is present with Udo = 10^. 



